For each i 



/ * rows * / 



for each j 
for each u 



/* columns */ 
/* u = 0,1, . . . , 



k */ 



MM u (i,j) = Mu(i,j)+ 2wu(x 9 y)e u (i-x 9 j-y) 
endfor (u) 

ExOut(i,j) = argmin (2 v u |MM u (k,l) -c u | p ) 1/p 

/* c u is the (u+l)-th coordinate of c 
(e 0 {i,j), ... e k (i,j)) = (MMo(i,j), MMi(i,j), . 



endfor (j) 
endfor (i) 

Set embedded source image M 0 ' as the first 
coordinates of ExOut. 



MM k (i, j) ) - ExOut (i, j) 



Figure 1A 
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For each iteration /* iteration */ 

for each i /* rows */ 

for each j /* column */ 

for each member d of C /* search through 

all possible members 
of C */ 

set ExOut (i, j) = d 
compute 

v(d) = |L(ExOut 0 - Mo) | 2 +|L(PExOut! - M 0 )| 2 
endfor (d) 

set ExOut (i,j) = argmin v(d) 

d 

endfor (j) 
endfor (i) 

endfor (iteration) or until ExOut has not changed 
between two consecutive iterations. 
Set embedded source image M 0 ' as the first 
coordinates of ExOut. 



Figure 5 A 
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